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Method, System, and Computer Program Product for Enhanced 
Resolution, Automatically- Calibrated Magnetic Position Sensor 



Background 

The present invention method, system, and computer program product provides, among other 
things, improved magnetic position sensor. Position sensors that use Hall effect elements to 
detect periodic magnetic discontinuities are well known and have proven useful in industrial 
applications related to aluminum die casting. Also, several techniques to improve the 
resolution of these sensors are known. However, existing forms of this sensor suffer from a 
combination of the following deficiencies, but not limited thereto: 

• Many require large chains of analog electronic circuitry, increasing expense and 
increasing the sensor's sensitivity to electromagnetic noise and environmental changes. 

• Many circuits neglect to adequately compensate for changes in amplitude from the Hall 
effect elements. 

• Many systems require manual calibration of analog components like potentiometers. 

• Many designs require the size and complexity of the circuitry to scale with the expected 
resolution enhancement. For example, a further doubling of resolution requires a 
doubling in size and complexity of circuit components. 

• Many designs require operations that are difficult and costly to implement, such as 
digital division operations. 



Brief Summary of the invention 

The present invention is a complete magnetic position sensor and related method and 
computer program code, which overcomes the problems suffered by traditional and improved 
magnetic positions sensors of this variety. It is capable of increasing resolution by an arbitrary 
degree, limited only by input signal noise. The sensor has proven capable of resolutions of 
about ,0015 inch, 32 times the resolution of the underlying magnetic stripes, operating at 500 
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inches per second. The s>stem is small, simple, efficient, and almost entirely digital. It 
requires minimal analog circuitry and calibrates itself automatically. 



Detailed Description of the Invention 

The exemplary individual parts of the invention are described below. 
Hydraulic Ram 

A hydraulic ram, with regularly spaced bands of magnetic material around the girth of the ram, 
is identified in Figure 1. These bands are spaced at intervals of 1/20* of an inch, although 
other intervals are possible. In the text that follows, this distance is indicated by the variable 
toothPeriod, 

Hall Effect Sensors 

Four Hall effect magnetic sensors, labeled Halll, Hall2, HalLJ, and Hall4 in the Figure 1, are 
mounted in close proximity to the ram. These sensors are spaced at precise distances with 
respect to one another, such that the first sensor is at 0 degrees, the second at 90 degrees, the 
third sensor at 180 degrees, and the fourth sensor at 270 degrees. These angular 
measurements are relative to the toothPerioi where 360 degrees is equal to one toothPeriod, as in 
the following equation: 

i position 
angle = 360 x- 



toothPeriod 

As the ram moves, the four hall sensors are fixed in position, such that the alternating bands of 
magnetic material pass beside the Hall effect sensors. Thus the Hall effect sensors provide 
four time- varying, sinusoidal Hall voltages, with precise phase relationships. If variable position 
is assigned to be the linear position of the ram (where the origin, pcsition=Q, is assigned to be 
directly overtop a tooth), then the following relationships exist between p and the hall voltages: 

hl = cosf 3 60°x P° sition ) 
\ toothPeriod J 



h2 = cosf360° X p0Siti ° n - -90-1 
^ toothPeriod J 
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h3 = co/360°x position -ISO') 
^ toothPeriod J 

h4 = cosf360'x P° sition -270°] 
V toothPeriod J 

Differential Amplification 

Hall voltages hi and h3 are passed through a differential amplifier (labeled AmpO in Figure 1), 
which provides both amplification and common- mode noise rejection. The output voltage for 
AmpO, labeled aO, follows the equation: 

aO = gain x (input + - input, ) 
Thus, the output can be expressed: 

'cosf 360- x P° sition 1-cosf 360° x P ° siti ° n -180°" 
toothPeriod J I toothPeriod 



aO = gain x 

Since cos(x + 1 80° ) = - cos(x) } f or all x: 

a0 = 2xgainxcosf360 o x P° sition ) 
\ toothPeriod ) 

Likewise, signals h2 and h4 are connected to Amp90, yielding: 

a9Q = 2x gain xsinf 360° x P Qsitl0n ) 
{ toothPeriod J 

Digital Controller 

The digital controller manages the flow-of-control for the components in the enclosing box in 
Figure 1. For simplicity, wires that connect the digital controller to the individual components 
have not been drawn. The invention uses a MSP430 Mixed Signal Processor manufactured by 
Texas Instruments, although many other choices are available, including Field Programmable 
Gate Arrays. 

Offset Calibration 

The amplified signals, aO and a90, have an unpredictable DC offset component due to 
manufacturing tolerances, placement tolerances, and variation in magnetic flux intensity. 
However, offsets must be corrected to a fixed level of 1.5V to assure that aO and a90 are 
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correctly sampled by the analog- to-digital circuitry, precise dc offset is attained using a simple 
control loop. For signal aO, the control loop consists of the analog-to-digital component 
ADOO, the averaging component AvgO, the digital- to- analog component DAQ), and the 
subtracter SubO. Calibration is performed only when power if first applied to the sensor, a 
process managed by the digital controller. During calibration, a 0 volt signal is initially 
outputted by AvgO. Thus, the output from SubO follows its input, aO. The ADGO then 
samples aO at 100 kilosamples per second, and AvgO performs a time- average of these samples 
over a .25 second interval. The component's output, labeled errO in Figure 1, is generated by 
subtracting the resultant average value from one 1.5 volts: 

err0 = 1.5 V -average 

, and outputted digitally by the AvgO component. At this point, the subtractor will remove this 
error component from subsequent analog values of aO. Hereafter, the calibrated signal sO will 
have a precise DC offset of 1.5 volts. By the same process, signal s90 will also have a correct 
DC offset of 1.5 volts. 

While AvgO and Avg90 could be implemented as separate external components or internal 
blocks of an FPGA, this invention implements them in software within the MSP430 digital 
controller. 

Measurement 

The invention contains two analog- to-digital converters, labeled ADOO and ADC90 in Figure 
1, which translate analog signals within the range of 0V to 3V to 12-bit signed integer values 
within the range -2048 to 2047. After conversion, calibrated signals dO and d90 vary 
sinusoidally with respect to the ram position. These signals are centered at a DC offset of 0, 
the middle of the 12-bit sampling window. 

While ADC0 and ADC90 could be implemented as separate external components, this 
invention uses the analog-to-digital converters that are available on the MSP430 digital 
controller. 

Signal Selection 

After measurement, digital signals dO and d90 are passed to signal selection component label 
Sel in Figure 1. This component simply compares the absolute magnitudes of dO and d90 and 
returns the signal with the lesser magnitude as its primary output, labeled p in the diagram, and 
returns the signal with greater magnitude as hs secondary output, labeled s in the diagram. The 
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output signal labeled sel in Figure 1 is a digital signal generated inside the component: it is zero 
when dO is the primary signal and one when d90 is the primary signal. 

While the signal selection component could be implemented as a separate external component, 
this invention implements it in software within the MSP430 digital controller. 

Normalization 

The normalization component estimates the amplitude of the incoming primary and secondary 
signals. These signals, p and s, have the form amplitude x cos(x) m d amplitudex sin(x) or 
vice- versa. Thus, according to the geometric identity sin 2 (x) + cos 2 (x) = 1 : 

^p 2 +s 2 = -^/amplitude 2 x (cos 2 (x) + sin 2 (x)) = amplitude 
Rewritten, the equation becomes 

amplitude = -^p 2 + s 2 

This appears to require a digital square- root operation. However, since the signals p and s are 
represented by twelve- bit integers, only 11-bit amplitudes are meaningful. Thus, the operation 
is gready simplified through the use of a lookup table. The lookup table works as follows: 
Elements are stored in a read-only memory with addresses 1,2,... 2048. The value at each 
memory address is equal to the square of the address. For example the value at address 4 is 16, 

the value at address 5 is 25. Tlius, for a given square x , the task of finding Vx is reduced to a 
11-step binary search. 

After the amplitude is estimated, the normalization component then scales the primary signal, 
p according to the following equation: 

2048 xp 

np = — 

amplitude 

, where np is the output of the component, labeled in Figure 1. Again, a lookup table is used 
to eliminate the costly division operation. The table works as follows: Elements are stored in a 
read-only memory with addresses 1,2,... 2048, each address corresponding to a possible 
amplitude as calculated above. The value at each memory address is equal to: 
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value = 



2048 . 



address 



The resulting memory value is then multiplied by p and the result is shifted right by 10 bits. 
This final 'right shift' operation is a fast, digital technique to divide by 2 10 . Thus, the equation 
for np becomes: 



np = 



2048 



amplitude 



-x2 



10 



O .,o 2048 xp 
xpx2 a — 

amplitude 



This technique maps signal p, a value along a sinusoidal signal with unknoun amplitude, into a 
value along a sinusoidal signal with knom amplitude, 2048, While the normalization 
component could be implemented as a separate external component, this invention 
implements it in software within the MSP430 digital controller. 

Angle Table 

The angle table is a read-only memory that maps sinusoidal input signal np into a phase angle, 
a discrete interval within the phase number line. Hie angle table follows this equation 



angle = 



resolution ( np ^ 

x arccos — — 

360° \204%) 



, where resolution is the number of subintervals within the range [0, 360°]. Constant resolution 
ultimately dictates the resolution of the sensor; it represents the number of discernable output 
events for every tootkPeriocL This invention uses a resolution of 32, although any other 
resolution constant would suffice. 

The equation listed above requires the arccosine operation, a function whose domain is [-1,1] 
and range is [0,180°]. The expression (np/2048) maps the signal np, which varies from[- 
2048,2047], to the necessary [1,1] interval. The multiplication by (resolution / 360°) constant 
and the subsequent floor operation maps the result of the arccosine operation onto the integer 
interval [0, (resolution/2) - 1]. 

Although the above equation describes the nature of the operations, the implementation is 
greatly simplified through the use of a lookup table, which eliminates the cosdy arccosine and 
division operations. The table works as follows: Elements are stored in a read-only memory 
with addresses -2048,-2047,... 2046, 2047 corresponding to possible values of the signal np. 
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The value at each memory address is equal to the result of the entire equation defined above. 
In other words, the entire equation is reduced to single memory retrieval operation. 

While the angle table could be implemented as a separate external component, this invention 
implements it in software within the MSP430 digital controller. 

Phase Translation 

The phase translation component performs two important corrections on signal ang: 

• Because the cosine function is not one-to-one, its inverse function arccosine only 
resolves to the range [0, (resolution/2) - 1], corresponding to the phase interval [0, 
180°]. Thus, a second step is necessary to correcdy resolve the phase angle to the 
interval [0, resolution -1] corresponding to the entire phase interval [0, 360°]. 

• If signal ang is derived from signal d90 (a sinusoid that is 90 degrees ahead of dO), the 
resulting phase angle must be shifted forward by 90 degrees. Figure 2 summarizes the 
relationships that exist between dO and d90. 

Both of these corrections are described below, the phase translation component's behavior is 
broken into several cases: 

Derived from dO 

If signal sd (an output from the signal selection component) is zero, then dO was selected as 
the primary signal and ang was derived from dO. As outline earlier, this case requires that 
signal ang be correcdy shifted into the interval [0°, 360°]. The corrected phase interval can be 
determined by inspecting the sign of signal d90, which appears at input s. If d90 is greater 

than 0, then the phase is correcdy situated in the [0, n] , as shown in Figure 3. Thus the 
output, phi, should match the input, ang: 

phi = ang 

If d90 is less than 0, then the phase should be corrected to be in the interval [rc,2n] j a$ shown 
in Figure 3. The output should be corrected according to the following equation: 

phi = resolution - 1 - ang 

As an example, consider Figure 3. The figure shows normalized sample of 900 appears at 
signal np. The angle table will map this to a phase angle of 64°, although it is possible that that 
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the correct angle is 296° (since there are two possible phase angles that correspond to value 
900). To determine the correct angle, the value of the secondary signal s is inspected. This is 
shown in Figure 4. If the value of $ is positive, then the correct angle is 64°. If it is negative, 
the correct angle is 360° - 64° = 296°. 

Derived from d90 

If sel equals one, then ang was derived from d90. Hie correct phase interval can be 
determined by inspecting the sign of signal dO, which appears at input s. If dO is less than zero, 
then only the 90 degree correction is necessary. This operation is described in the following 
equation 



phl = 



resolution > 

ang + 



% resolution 



4 

, where the % symbol represents modulo division. If dO is greater than zero, then two 

corrections are necessary. Signal ang is first mapped into the interval IX2ti] m Additionally, the 
signal is shifted forward by 90 degree correction: 

phi = ^(resolution - 1 - ang)+ reSQ ^ on % resolution 

After these corrections, signal phi correcdy represents the phase angle encoded as an unsigned 
integer between the range [0, resolution - 1]. While the phase translation component could be 
implemented as a separate external component, this invention implements it in software within 
the MSP430 digital controller. 

Phase Register 

Following the output phase, described later, that maris the end of each sensing iteration, the 
output of the phase register (ph2) is set to each the output of the phase translation component 
(phi). Thus, the phase register always holds the phase value from the previous iteration. 
While the phase register component could be implemented as a separate external component, 
this invention implements it in software within the MSP430 digital controller. 

Phase Subtractor 

This phase subtractor component determines the amount of angular movement for the 
current iteration, relative to the phase angle from the previous iteration. This movement is 
simply the difference between the new phase value and the phase value for the previous 
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iteration, as stored in the phase register. The resulting digital output from this component is 
the signal labeled mm in Figure 1. While the phase subtractor component could be 
implemented as a separate external component, this invention implements it in software within 
the MSP430 digital controller. 

Overflow Corrector 

The overflow corrector handles overflow and underflow cases. These situations arise when 
the phase angle transitions from 359° to 0° or from 0° to 359°. Such situations yield an 
apparent movement of -359° or 359°, respectively, though an actual movement of only 1° or - 
1°, respectively. The following technique corrects this erroneous movement. If the apparent 
movement, labeled win Figure 1, is greater than (resolution / 2), then m is corrected by the 
equation: 

mc = m- resolution 

If the apparent movement is less than (resolution /2), then m is corrected by the equation: 

mc = m + resolution 
If neither of the above conditions is true, then the input is unmodified: 

mc = m 

While the overflow corrector could be implemented as a separate external component, this 
invention implements it in software within the MSP430 digital controller. 

Output 

The output component generates three output signals labeled lead, trails and vdidm Figure 1. 
These signals are toggled at the end of each sensing iteration. First, the lead and trail signals are 
toggled to emulate a traditional quadrature output (which is usually generated by zero crossing 
detection circuits from sinusoids in quadrature). In this invention, a simple finite state 
machine controls the lead and trail signals, as described in Figure 5. The signal rrc determines 
the number of transitions to follow along the state machine, at 100 nanoseconds per state. If 
ttris positive, then "forward" transitions are followed; if ^ris negative, then "backward" 
transitions are followed. 

After this signaling is complete, the vdid signal is raised high for 100 nanoseconds, then 
returned back to zero. This signal is used to indicate to decoding logic that the output 
signaling has finished. While the output component could be implemented as a separate 
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external component, this invention implements it in software within the MSP430 digital 
controller. 

Still other embodiments will become readily apparent to those skilled in this art from 
reading the above- recited detailed description and drawings of certain exemplary 
embodiments. It should be understood that numerous variations, modifications, and 
additional embodiments are possible, and accordingly, all such variations, modifications, sizes, 
levels, and embodiments are to be regarded as being within the spirit and scope of this 
document. For example, regardless of the content of any portion (e.g., tide, section, abstract, 
drawing figure, etc.) of this application, unless clearly specified to the contrary, there is no 
requirement for any particular described or illustrated activity or element, any particular 
sequence of such activities, or any particular interrelationship of such elements. Moreover, any 
activity can be repeated, any activity can be performed by multiple entities, and/ or any element 
can be duplicated. Further, any activity or element can be excluded, the sequence of activities 
can vary, and/ or the interrelationship of elements can vary. Accordingly, the descriptions and 
drawings are to be regarded as illustrative in nature, and not as restrictive. 

The invention may be embodied in other specific forms without departing from the 
spirit or essential characteristics thereof. The foregoing embodiments are therefore to be 
considered in all respects illustrative rather than limiting of the invention described herein. . 

REFERENCES 

The following US. Patents are hereby incorporated by reference herein in their entirety. 

> 5012239, 

> 4630928, 

> 5442313, 

> 5067089, 

> 5719789, 

> 4587485, 
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> 4972080 

> 3,956,973 

> 5,041784, 



In summary, the present invention method, system, and computer program product 
provides, among other things, a magnetic position sensor for resolution-critical, harsh- 
environment industrial applications like aluminum die casting. It outperforms existing sensors 
in every category- it operates at resolutions up to about 32 times higher than existing 
sensors at speeds of up to about 500 inches per second. It is inexpensive and can 
automatically calibrates itself. 

The invention is useful for any conceivable position sensing application. For example, it is useful 
for sensing displacement along pneumatic or hydraulic cylinders. These applications are common in 
aluminum die casting, medical robotics, forestry and wood manufacturing, etc. 

The design of the present invention could be applied to any sensor that provides analog, 
quadrature signals. Such signals are found in any sensors, most notably optical encoders, etc. 
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Appendix: Source Code 

The following is C source code, used to program the MSP430F149 mixed signal processor. The source code 
consists of n files - main.c, cosine.c, and cosineconstants.c. 

main.c 

#include <msp430xl4x.h> 
tdefine COSINERES 32 

fdefine F \ 

switch (quadrature) { \ 

case 0: quadrature = 1; break; \ 

case 1: quadrature = 3; break; \ 

case 2: quadrature = 0; break; \ 

case 3: quadrature = 2; break; \ 
} \ 

OUT = quadrature; \ 
offset++; 

#define B \ 

switch (quadrature) { \ 

case 0: quadrature = 2; break; \ 

case 1: quadrature = 0; break; \ 

case 2: quadrature = 3; break; \ 

case 3: quadrature = 1; break; \ 
} \ 

P40UT = quadrature; \ 
offset--; 

void algorithm (void) ; 
long offset; 
int zeroA; 
int zeroB; 

unsigned char quadrature = 0; 

void calibrate (void) { 
long i; 

long aveA = Ob- 
long aveB = Ob- 
long iter = 500000; 
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LCDDISPLAY ("Calibrating. . .", 14) ; 

LCDSECONDLINEO ; 

LCDDI SPLAY ("Move Sensors", 12); 



for (i=0; Kiter; i++) { 
aveA += ADC12MEM0; 
aveB += ADC12MEM1; 

} 

zeroA = aveA / iter; 
zeroB = aveB / iter; 

LCDCLEAR () ; 



void main (void) { 

// temporary variable used during calibrations 
int calib; 

// only use the watchdog timer if we're using the LCD 
lifdef LCD 

WDTCTL = WDT_ADLY_1000; // WDT ls/4 interval timer 

IE1 |= WDTIE; // Enable WDT interrupt 

#else 

WDTCTL = WDTPW + WDTHOLD; // Stop watchdog timer 

lendif 

// setup the clock for fast as she'll go 
DCOCTL = DCO0 + DCOl + DC02; 
BCSCTL1 |= RSELO + RSEL1 + RSEL2; 

LCDINITIALIZEO ; 



// PWM 

TACTL = TASSEL1 + TACLR; 
CCRO = 4095; 
CCTL1 = OUTMOD_7; 
CCR1 = 2048; 
CCTL2 = OUTMOD 7; 



// SMCLK, Clear Tar 

// PWM Period 

// CCR1 reset/set 

// CCR1 PWM duty cycle 

// CCR2 reset/set 
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CCR2 = 2048; 
P1DIR |= OxOC; 
P1SEL |= OxOC; 
TACTL |= MCO; 



// CCR2 PWM duty cycle 

// PI. 2 and PI . 3 output 

// PI. 2 and PI. 3 TA1/2 otions 

// Start Timer_A in up mode 



// ADC 
// 



P6SEL |= 0x88; 

ADC12CTL0 = ADC12ON+SHT0_2+MSC; 
ADC12CTL1 = SHP + ADC12SSEL_2; 
ADC12CTL1 |= CONSEQ_3; 
ADC12MCTL0 * INCH_7; 
ADC12MCTL1 - INCH_3 + EOS; 
ADC12CTL0 |= ENC; 
ADC12CTL0 |= ADC12SC; 



// Enable A/D channel AO 

// Turn on ADC12, set sampling time 

// Use MCLK for conversion timer. 

// Use consecutive mode 



// Enable conversions 
// Start conversion 



// Port 

// 

P4DIR | = OxfF; 

// Go 

// 

calibrate () ; 

calib = CCR1 - (zeroA - 2048); 
if (calib < 0) calib = 0; 
if <calib > 4095) calib = 4095; 
CCR1 = calib; 

calib = CCR2 - (zeroB - 2048); 
if (calib < 0) calib = 0; 
if (calib > 4095) calib = 4095; 
CCR2 = calib; 

algorithm () ; 

> 



14 



00949-01 



cosine.c 

tinclude "main.h" 

#define READ \ 

a = ((int) ADC 12 MEMO) - zeroA;\ 

b = ((int) ADC12MEM1) - zeroB; \ 

if (a < 0) { aMag = a * -1; } else { aMag = a; ) \ 

if (b < 0) { bMag = b * -1; } else { bMag = b; ) 

extern const long div2 [] ; 
extern const long squares []; 
extern const unsigned char cosine []; 
extern unsigned char quadratures- 
extern long offset; 
extern int zeroA; 
extern int zeroB; 
long a; 
long b; 

int aMag, bMag; 

void algorithm (void) { 
signed char last = 0; 
signed char delta = 0; 
signed char current = 0; 
int arrayPos; 
long amp = 800; 
long ampSQ; 
long outstanding = 0; 

int 1; 
int r; 

while (1) { 

P40UT A = 0x10; 
READ; 

ampSQ = a * a; 
ampSQ += b * b; 

// This next block of code performs a cheap square root 
if (ampSQ < squares [1024 ] ) { 
1 = 0; r = 1023; 
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} else { 

1 = 1024; r = 2047; 

} 

if (ampSQ < squares [1+512] ) 

r -= 512; 
} else { 
1 += 512; 

} 

if (ampSQ < squares [1+256] ) 

r -= 256; 
} else { 
1 += 256; 

} 

if (ampSQ < squares [ 1+128] ) 

r -= 128; 
} else { 
1 += 128; 

} 

if (ampSQ < squares [ 1+64 ] ) 

r -= 64; 
) else ( 
1 += 64; 

) 

if (ampSQ < squares [ 1+32] ) 

r -= 32; 
} else ( 
1 += 32; 

} 

if (ampSQ < squares [1+16] ) 

r -= 16; 
} else { 
1 += 16; 

) 

if (ampSQ < squares [1+8 ] ) { 

r -= 8; 
} else { 
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1 += 8; 

} 

if (ampSQ < squares [1+4 ] ) { 

r -= 4; 
} else { 
1 += 4; 

} 

if (ampSQ < squares [ 1+2 ] ) { 

r -= 2; 
) else ( 
1 += 2; 

} 

if (ampSQ — squares [1]) { 

amp = 1; 
} else { 

amp = 1+1; 

} 

// We use either signal a or signal b as the "lookup" source, 
// depending on which signal is smaller (smaller is better) 
if (aMag < bMag) { 

arrayPos = ( (a * div2[amp]) » 10) + 2048; 

if (arrayPos < 0) arrayPos = 0; 

if (arrayPos > 4095) arrayPos = 4095; 

current - cosine [arrayPos] ; 

if (b < 0) current = COSINERES - 1 - current; 
} else { 

arrayPos = ( (b * div2[amp]) » 10) + 2048; 

if (arrayPos < 0) arrayPos = 0; 

if (arrayPos > 4095) arrayPos = 4095; 

current = cosine [arrayPos] ; 

if (a > 0) current = COSINERES - 1 - current; 
current = (current + (COSINERES / 4)) % COSINERES; 

» 

delta = current - last; 
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// check for overflow or underflow 
if (delta > (COSINERES / 2)) { 

delta -= COSINERES; 
} else ( 

if (delta < - (COSINERES / 2) ) { 
delta += COSINERES; 

} 

} 

switch (delta) { 



case 


-16 


: B; B; B; B; B; B; B; B; B; B; 


B; 


B; B; B; B; B; break 


case 


-15 


: B; B; B; B; B; B; B; B; B; B; 


B; 


B; B; B; B; break; 


case 


-14 


: B; B; B; B; B; B; B; B; B; B; 


B; 


B; B; B; break; 


case 


-13 


: B; B; B; B; B; B ; B; B; B; B; 


B; 


B; B; break; 


case 


-12 


: B; B; B; B; B; B; B; B; B; 3; 


B; 


B; break; 


case 


-11 


: B; B; B; B; B; B; B; B; B; B; 


B; 


break; 


case 


-10 


: B; B; B; B ; B; B; 3 ; B; B; B; 


break; 


case 


-9: 


B; 


B; 


B; B; B; B; B; B; B; break; 




case 


-8: 


B; 


B; 


B; B; B; B; B; B; break; 






case 


-7: 


B; 


B; 


B; B; B; B; B; break; 






case 


-6: 


B; 


B; 


B; B; B; B; break; 






case 


-5: 


B; 


B; 


B; B; B; break; 






case 


-4 : 


B; 


B; 


B; B; break; 






case 


_ o . 

"j! 


B; 


B; 


B; break; 






case 


-2: 


B; 


B; 


break; 






case 


-1: 


B; 


break; 






case 


1: 


F; 


break; 






case 


2: 


F; 


F; 


break; 






case 


3: 


F; 


F; 


F; break; 






case 


4: 


F; 


F; 


F; F; break; 






case 


5: 


F; 


F; 


F; F; F; break; 






case 


6: 


F; 


F; 


F; F; F; F; break; 






case 


7: 


F; 


F; 


F; F; F; F; F; break; 






case 


8: 


F; 


F; 


F; F; F; F; F; F; break; 






case 


9: 


F; 


F; 


F; F; F; F; F; F; F; break; 




case 


10: 


F; 


F; 


F; F; F; F; F; F; F; F; 


break; 


case 


11: 


F; 


F; 


F; F; F; F; F; F; F; F; 


F; 


break; 


case 


12: 


F; 


F; 


F; F; F; F; F; F; F; F; 


F; 


F; break; 


case 


13: 


F; 


F; 


F; F; F; F; F; F; F; F; 


F; 


F; F; break; 


case 


14: 


F; 


F; 


F; F; F; F; F; F; F; F; 


F; 


F; F; F; break; 


case 


15: 


F; 


F; 


F; F; F; F; F; F; F; F; 


F; 


F; F; F; F; break; 


case 


16: 


F; 


F; 


F; F; F; F; F; F; F; F; 


F; 


F; F; F; F; F; break; 
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} 

// raise, then lower the valid line 

P40UT |= 0x04; 

_NOP(); 

P40UT S= ! (0x04) ; 
last = current; 

} 

) 



00949-01 



Cosine Constants 

Unlike the previous code written in the C language, the following code is written in Java. Once executed, it 
writes a large collection of tables to the screen or standard output. This resulting text is then saved in a C file 
called cosineconstants.c. These are the tables necessary to compile cosine.c, listed above. 

public static void main (String [ ] args) { 
String s = ("unsigned char cosine [] = {"); 
for (double i=0; i < 4096; i++) { 

double y = (i - 2048d) / 2048d; // returns 0 to pi 
double x = Math. floor ( (Math . acos (y) / Math. PI) * 16d )/ 
s+= ((int) x) + ", »; 

// if this is a newline 

if (<(i + 1) % 64 == 0) & (i ! = 0)) 

s += "\n n ; 

} 

System. out. println(s) ; 



s = ("float divisorst) = {W); 

for (float i=0; i < 2048; i++) { 

String c = + (int) (<2048f / i) * 1024); 

for (int j = c. length (); j < 10; j++) { c += " ) 

s+= c + ", "; 

// if this is a newline 

if (((i + 1) % 64 == 0) & (i != 0)) 

s += "\n"; 

} 

System. out. println (s) ; 

s = ("const long squares [] = (\n M ); 
for (long i=0; i < 2048; i++) { 
String c = ,,n + i*i; 

for (int j = c.length(); j < 10; j++) { c += " "; } 
s+= c + ", "; 

//if this is a newline 

if ({(i + 1) % 64 == 0) & (i != 0)) 

s += "\n"; 

} 

System. out .println (s) ; 
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